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AIRFOIL  PROFILES  FOR  MINIMUM  PRESSURE  DRAG  AT  SUPERSONIC 
VELOCITIES  -  APPLICATION  OF  SHOCK-EXPANSION  THEORY, 
INCLUDING- CONSIDERATION  OF  HYPERSONIC  RANGE 
By  Dean  R.  Chapman 


SUMMARY 


A  theoretical  investigation  is  made  of  airfoil  profiles  at  supersonic 
velocities  to  determine  the  shapes  having  minimum  pressure  drag  at  zero 
lift  for  various  given  auxiliary  conditions.  Shock-expansion  theory  is 
employed,  thereby  extending  the  applicability  of  the  results  through  the 
hypersonic  range.  Curves  are  presented  for  Mach  numbers  of  1.5,  2,  3, 
k,  6,  8,  and  00  which  enable  the  shape  and  the  drag  of  an  optimum  profile 
to  be  determined  readily  if  the  base  pressure  is  known  from  experiments. 
Examples  are  presented  of  optimum  profiles  determined  with  the  aid  of 
experimental  base  pressure  data.  Variations  in  profile  shape  are  inves¬ 
tigated  to  provide  information  on  the  degree  to  which  deviations  in 
shape  from  the  optimum  can  be  made  without  resulting  in  a  significant 
drag  increase. 

A  comparison  of  optimum  profiles  determined  by  the  shock-expansion 
method  of  this  report  with  corresponding  profiles  determined  by  the 
linearized -theory  method  of  a  previous  report  shows  only  small  differ¬ 
ences  in  shape  at  Mach  numbers  up  to  infinity  even  though  the  linearized 
theory  at  high  supersonic  Mach  numbers  breaks  down  completely  insofar  as 
the  drag  of  the  profile  is  concerned.  The  experimentally  observed 
dependence  of  base  pressure  on  trailing-edge  thickness  is  found  to  have 
a  significant  effect  on  the  shape  and  drag  of  optimum  profiles  of  small 
thickness  ratio.  Curves  are  presented  which  show  that  for  thin  airfoils 
the  use  of  a  trailing-edge  thickness  considerably  greater  than  the  theo¬ 
retical  optimum  can  result  in  an  excessive  drag  penalty  at  moderate 
supersonic  Mach  numbers,  though  not  at  hypersonic  Mach  numbers. 


INTRODUCTION 

In  1933  Saenger  observed  that  for  the  extreme  case  of  flow  at 
infinite  Mach  number  an  airfoil  designed  to  have  minimum  pressure  drag 
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would  have  its  maximum  thickness  at  the  trailing  edge.  (See  refer¬ 
ence  1.)  A  related  result  can  be  inferred  from  the  numerical  calcula¬ 
tions  of  Ivey  (reference  2)  which  indicate  that  the  drag  of  a  10-percent- 
thick-wedge  airfoil  at  a  Mach  number  of  8  is  less  than  that  of  a  double¬ 
wedge  airfoil  having  the  same  thickness  ratio.  In  both  of  these  cases 
the  desirability  of  employing  a  thick  trailing  edge  in  conjunction  with 
a  small  surface  slope  may  be  attributed  to  the  fact  that  at  hypersonic 
Mach  numbers  the  suction  forces  (forces  due  to  pressures  below  ambient) 
are  small  compared  to  the  positive  pressure  forces,  even  when  the  suction 
force  corresponds  to  a  vacuum.  Recently,  Smelt  (reference  3)  has  dis¬ 
cussed  this  latter  characteristic  of  hypersonic  flow  and  its  possible 
application  to  the  determination  of  efficient  airfoil  shapes  for  use  at 
very  high  Mach  numbers.  The  investigations  of  Saenger,  Ivey,  and  Smelt, 
however,  do  not  provide  general  quantitative  information  on  the  airfoil 
profile  having  minimum  pressure  drag  in  hypersonic  flow  because  of  two 
limitations  of  their  analyses:  Airfoils  having  a  trailing-edge  thickness 
less  than  the  maximum  airfoil  thickness  were  not  considered,  and  the 
airfoil  structural  characteristics  were  not  considered  (comparisons  were 
made  on  the  basis  of  a  given  airfoil  thickness). 

At  low  and  moderate  supersonic  Mach  numbers  the  suction  forces  on 
an  airfoil  can  amount  to  several  times  the  positive  pressure  forces, 
particularly  if  an  excessively  thick  trailing  edge  is  employed.  As  a 
result,  the  optimum  trailing-edge  thickness  in  this  lower  Mach  number 
range  depends  to  a  great  extent  on  the  base  pressure.  By  presuming  that 
the  base  pressure  is  known  from  experiments,  and  that  the  airfoil  profile 
must  satisfy  a  given  structural  requirement  (such  as  a  given  section 
modulus  or  a  given  section  moment  of  inertia),  a  method  of  calculating 
the  profile  of  minimum  pressure  drag  at  zero  lift  has  been  developed  by 
the  present  writer  in  reference  4.  Although  the  basic  equations  devel¬ 
oped  in  reference  4  for  calculating  such  profiles  are  applicable  to 
higher-order  theories,  a  detailed  solution  was  given  only  for  the  case 
of  linearized  supersonic  flow. 

Because  of  the  well-known  shortcomings  of  linearized  theory,  it  was 
thought  worthwhile  to  conduct  an  investigation  parallel  to  that  of 
reference  4  in  which  the  shape  and  drag  of  optimum  airfoils  are  determined 
from  shock-expansion  theory  instead  of  linearized  theory.  Shock-expansion 
theory  appears  adequate  for  this  purpose,  particularly  in  view  of  the 
recent  investigation  of  Eggers  and  Syvertson  (reference  5)  which  indicates 
that  shock-expansion  calculations  accurately  determine  surface  pressures 
on  thin  airfoils  in  inviscid  flow  at  Mach  numbers  from  just  above  that 
for  bow-wave  attachment  to  infinity.  From  the  viewpoint  of  the  engineer 
who  always  has  to  make  design  compromises,  it  was  thought  desirable  in 
the  present  study  also  to  determine  how  much  the  optimum-profile  shape 
can  be  altered,  especially  near  the  trailing  edge,  and  still  not  increase 
the  drag  excessively.  The  purposes  of  the  present  investigation,  there¬ 
fore,  were  (l)  to  develop  a  usable  method  for  determining  the  shape  and 
drag  of  optimum  profiles  in  the  Mach  number  range  beyond  that  covered 
adequately  by  linearized  theory  (step-by-step  details  involved  in  applying 
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the  method  developed  are  given  in  an  appendix),  and  (2)  to  determine 
curves  showing  the  rate  at  which  the  total  pressure  drag  increases  as 
the  profile  shape  deviates  from  the  optimum. 


NOMENCLATURE 


cd 

C 

h 

H 

I 

T 

Jn 

l 

L 

M 

n 

P 

Pt 

P 

Re 


t 

V 


airfoil  chord 
pressure  drag  coefficient 
constant  depending  on  y 
trailing-edge  thickness 


(t/2) 


n-u 


dx 


dimensionless  trailing-edge  thickness  ^  J 

given  value  of  auxiliary  integral  c 

I 

dimensionless  value  of  I 


( l 


normalizing  factor  for  |(Y),  defined  by  equation  (10) 
length  of  surface  of  constant  thickness 
dimensionless  length  of  surface  of  constant  thickness 
Mach  number 

parameter  appearing  in  definition  of  I 
pressure 
total  pressure 

/  P  “  P  oo\ 

pressure  coefficient  /  ■j - 


gPccV* 


Reynolds  number 


chordwise  distance  from  leading  edge  to  first  downstream  position 
of  maximum  thickness 

maximum  thickness  of  airfoil 

velocity 

chordwise  distance  from  leading  edge  to  point  on  airfoil  surface 
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h 

y  ordinate  of  upper  surface  of  airfoil 

X  dimensionless  distance 

Y  dimensionless  distance 

y  ratio  of  specific  heats  (1.40  for  air) 

\  arbitrary  constant 

p  mass  density 

S  local  angle  of  inclination  of  airfoil  surface  with  respect  to 
chord  line  (tan“1y' ) 

£(Y)  characteristic  function  determining  optimum-profile  shape,  defined 
by  equation  (9) 

q(Y)  characteristic  function  \defined  by  equation  (12) 

Subscripts 

o  airfoil  surface  at  leading  edge 

1  airfoil  surface  at  trailing  edge 

00  free  stream 

b  base,  or  trailing  edge,  01  airfoil 

ca  circular-arc  biconvex  airfoil  having  sharp  trailing  edge 

Superscripts 

*  differentiation  with  respect  to  x 

ANALYSIS 

Solution  for  Arbitrary  Structural  Requirement  v 

As  in  reference  k,  it  is  assumed  throughout  this  analysis  that  the 
optimum  airfoil  has  a  sharp  leading  edge,  a  fixed  chord  length,  and  is 
set  at  the  zero-lift  angle.  It  also  is  assumed  that  the  flow  is  a  . 
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purely  supersonic  two-dimensional  flow  of  an  inviscid,  nonconducting, 
perfect  gas.  Since  the  surface  pressures  on  the  top  and  bottom  of  an 
airfoil  can  be  calculated  independently  in  a  supersonic  flow,  it  follows 
that  at  zero  lift  the  optimum  profile  will  be  symmetrical  about  the 
chord  line.  The  mathematical  problem  formulated  is  to  find  the  airfoil 
ordinate  function  y(x)  which  minimizes  the  pressure  drag, 

2  ,c  h  ,  . 

cd  =  c  4  Py’  dx  “  Pb  c  C1) 

for  a  given  value  of  the  auxiliary  integral 


I 


_  1  fc  y11 

C  J°  (t/2)a 


dx 


(2) 


By  selecting  various  values  of  the  parameters  n  and  cr,  a  wide  variety 
of  structural  requirements  can  be  represented  for  both  thin-skin  and 
solid-section  structures.  Some  of  the  different  structural  criteria 
represented  by  equation  (2)  are: 


n 

a 

Structural  criteria 

1 

0 

given  torsional  stiffness,  or  torsional  strength,  of 
thin-skin  structure  (given  cross-section  area) 

2 

0 

given  bending  stiffness  of  thin-skin  structure 

3 

0 

given  bending  stiffness  or  given  torsional  stiffness  of 
solid-section  structure 

2 

1 

given  bending  strength  of  thin-skin  structure 

3 

1 

given  bending  strength  of  solid-section  structure 

Basic  equations.-  The  equations  which  the  optimum-airfoil  ordinate 
function  y(x)  must  satisfy  can  be  obtained  by  considering  an  infinitesi¬ 
mal  variation  in  profile  shape  8y(x)  that  is  arbitrary  except  for  the 
requirement  that  81  =  0.  By  also  requiring  that  8c<j  =  0  the  following 
three  equations  result  (see  reference  4): 


d_ 

dx 


+  X,ny 


n-  1 


=  0 


(3) 


Pb  =  Pi 


w 


—  =  —  I  (I  =  l/(t/2)n-a) 
c  n 


In  the  derivation  of  these  equations  it  was  assumed  that  P  =  P(y',  yoS  Mj, 
but  no  particular  functional  form  was  assumed.  Equation  (3)  is  the  differ¬ 
ential  equation  which  the  airfoil  ordinate  function  y(x)  must  satisfy  along 
the  curved  surfaces,  for  example,  0A  and  BC  in  figure  1.  A  first  integral 
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of  this  differential  equation,  satisfying  the  condition  y*  =  0  at 
y  =  t/2,  can  be  obtained  by  multiplying  both  sides  by  y’  and  integrating. 

J*  SF  =  X  [  (t/2)n  -  7n]  (3a) 

Equations  (3)  and  (3a)  involving  the  arbitrary  constant  \  do  not  apply 
to  the  straight  midsection  (AB  in  fig.  l)  along  which  the  airfoil  thick¬ 
ness  is  constant.  Equation  (4),  termed  an  end  condition,  represents  the 
relation  which  must  be  satisfied  between  the  base  pressure  coeffi¬ 
cient  Pb,  the  surface  pressure  coefficient  P-j^  just  upstream  of  the 
trailing  edge,  and  the  corresponding  surface  slope  yx* .  As  will  be 
seen  this  equation  determines  the  optimum  trailing-edge  thickness. 
Equation  (5)  relates  the  optimum  length  of  straight  midsection  1  to 
the  dimensionless  value  of  the  structural  integral  I.  This  latter 
equation  shows  that  the  length  of  straight  mid section  is  always  zero 
when  the  auxiliary  condition  represents  a  given  stiffness  (a  =  0),  but 
for  the  values  of  n  considered  is  a  sizable  fraction  of  the  chord 
when  the  auxiliary  condition  represents  a  given  strength  (a  =  l). 


Solution  for  Y,  H,  and  X.-  From  equation  (3a)  the  constant  \  is 
readily  evaluated  in  terms  of  y0*.  There  results 


\(t/2)n  -  yQ»  {^yt)  (3b) 

By  employing  the  dimensionless  variables  Y  =  y/(t/2)  and  H  =  h/t, 
equations  (3a)  and  (3b)  yield 


yn  =  q_  y'^dp/dy* ) 

“  y0,2(Sp/8y')o 


(6) 


and  hence 

tjH  _  1  _  yi,a(dP/dy* )i 

"  y0,2(dp/dy')o 

It  is  to  be  noted  that  the  structure  of  these  equations,  and  all  sub¬ 
sequent  ones,  is  such  that  the  quantities  Pt>,  I,  and  which  are  pre¬ 
sumed  to  be  given  do  not  appeal1  as  independent  variables.  Instead,  they 
are  related  parametrically  to  the  shape  and  drag  of  the  optimum  profile 
through  the  parameter  y* .  For  example,  equation  (7)  gives 
HCYxS  Jo'  >  Mj,  and  equation  h  gives  Pt,(yi',  Jo'>  Mra),  hence  by  select¬ 
ing  arbitrary  values  of  yi '  the  function  H(Pb,  yo%  Mj  can  be  deter¬ 
mined.  Moreover,  as  will  be  seen  later,  I  depends  on  y0’ ,  M*,,  and  H, 
so  that  the  parametric  structure  of  the  equations  ultimately  yields  the 
desired  function  H(Pb,  I,  Mj. 


The  differential  equation  (3a)  and  the  appropriate  boundary  condi¬ 
tions  y( 0)  =  0,  y(s)  =  t/2,  and  y(s  +  l)  =  t/2  determine  x  as  a 
function  of  y  by  a  single  quadrature. 


<* 


* 
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J dp/8y'  dy 

./(t/2)'n  -  yE 


fep/dy«  dy 


n  _  vn 


(s  +  l)  + 


upstream 

surface 


8P/8y*  dy  on 

. downstream 

(t/2)n  -  y11  surface 


The  algebraic  sign  to  be  selected  for  each  radical  is  determined  by  the 
requirement  that  x  increases  as  y  increases  on  the  upstream  surface, 
and  as  y  decreases  on  the  downstream  surface.  The  above  relation  is 
more  conveniently  expressed  in  terms  of  the  dimensionless  quantities 
Y  =y/(t/2),  X  =  x/s,  and  L  =  l/ s.  After  eliminating  */~\  by  evaluating 


the  above  integral  at 
shape 


x  =  s,  there  is  obtained  for  the  .dimensionless 


i(Y)  ■  on  upstream  surface 

X  = 

L  +  |(Y)  on  downstream  surface 
where  £(Y)  is  defined  as  the  quotient  of  two  integrals1 

L  J  dp/dy'  Sl/J  1  -  Yh 


(8) 


l(Y)  =  ^- 

r 

Jo 


The  definite  integral  which  normalizes  |(y),  such  that 
for  sake  of  brevity  be  designated  by  jn.  Thus 


I (1)  =  1,  will 


Jn  = 


1  -  Yn 


Actually,  for  completeness  the  function  £(Y)  should  be  written  as’ 

|(Y;  n,  y0*)  because  it  depends  on  the  three  quantities  M^,  n,  and 
y0’  as  well  as  the  variable  Y.  For  brevity,  though,  it  is  written 
simply  as  £(Y). 

Since  x  is  a  double-valued  function  of  y  over  the  chord  length, 
£(Y)  is  also  a  double-valued  function  of  Y.  A  sketch  of  a  typical 
curve  of  Y  versus  |  is  shown  in  figure  2.  For  a  given  dimensionless 
ordinate  Y,  one  of  the  two  values  of  I  represents  the  dimensionless 
chordwise  distance  from  the  leading  edge  to  a  point  on  the  upstream 
surface,  whereas  the  other  value  represents  (apart  from  an  additive  % 
term  L)  the  dimensionless  chordwise  distance  from  the  leading  edge  to 
the  point  on  the  downstream  surface  which  has  the  same  ordinate  Y.  It 
is  to  be  understood  that  in  determining  X  from  equation  (8)  the  appro¬ 
priate  value  of  £(Y)  must  be  used  for  each  surface;  thus,  for  a  given  Y 

1Since  the  integrands  are  singular  at  the  point  Y  =  1,  numerical  compu¬ 
tations  of  the  function  £(y)  must  allow  for  this  singularity.  A 
simple  method  of  doing  this  is  outlined  in  appendix  A. 
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the  value  of  5 (Y)  appropriate  to  the  downstream  surface  always  is 
greater  than  for  the  upstream  surface. 

From  the  foregoing  and  the  fact  that  |(y)  is  independent  of  Pb 
and  or,  it  follows  that  a  curve  of  Y  versus  |(y),  such  as  is  illus¬ 
trated  in  figure  2,  determines  the  curved  portions  of  an  infinite  number 
of  optimum  profiles,  all  having  the  same  values  of  y0’,  M*,  and  n.  It 
maybe  noted  that  a  (equation  (5))  essentially  determines  the  length 
of  straight  midsection  which  is  to  be  placed  between  the  two  curved 
portions  after  separating  them  at  the  point  where  Y  =  1  (point  AB  in 
fig.  2);  whereas  Pb  (equations  (4)  and  (7))  essentially  determines  H, 
the  value  of  Y  beyond  which  the  downstream  portion  of  the  curve  is  not 
used  in  a  given  case.  It  is  noted  that  although  the  chord  is  fixed, 
the  value  of  £  corresponding  to  the  trailing  edge  is  not.  This  is 
because  |  =  x/s  changes  whenever  s  changes. 


As  may  be  deduced  from  equation  (7),  Yi*  determines  H  for  a 
given  n,  M^,  and  y0*.  Moreover,  y^  determines  Pb  from  equation  (4). 
Hence,  for  any  given  value  of  base  pressure  the  point  on  the  downstream 
surface  which  corresponds  to  the  trailing-edge  position  can  be  indicated 
on  each  £(Y)  curve.  (See  fig.  2  where  the  point  corresponding  to  zero 
base  pressure  is  indicated.) 


Solution  for  L,  I,  t/c,  and  s/c.-  Turning  now  to  the  determination 
of  the  optimum  length  of  straight  midsection  l,  one  sees  from  equa¬ 
tion  (5)  that  such  a  determination  will  also  give  I.  Since  Z  is  a 
function  of  y0’  and  the  given  quantities  n,  a,  Pb',  and  M^,  this 
enables  xy0’,  the  quantity  used  as  a  parameter  in  the  present  analysis, 
to  be  related  to  the  quantity  I,  which  is  a  more  convenient  one  to  use 
if  I  is  the  actual  quantity  given.  Starting  with  the  definition  of  T 
the  following  equations  result: 


I  =  l/(t/2)n"a 

=  - /G/Srn  dx  =  -J £/S  [1  -  (1  -  Yn)]  dx 

c  Jo  CJ0  ' 

-/H  (1  -  Yn)  g  dY 
c[_s  Jo  dY 

or,  by  using  the  relation 


dX  _  _1_  /dP/dy' 
dY  jn  J  1  -  Yn 


(ID 


which  follows  from  equations  (8),  (9),  and  (10),  there  results 


1  -  I 
s/c 


-  Yn) 


■  * 


U 


% 


jj* 


dY 
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For  convenience  the  right-hand  member  of  this  equation-  is  defined 
as  T)(H). 

n(H)  e  i  (1  -  Yn)  dY  (12) 

Jn  J0  *  oy' 

From  equations  (8)  and  (5), 


7  =  L  +  6(H) 


-§fT+  1(h) 


Combining  this  with  the  above  relation  between  I  and  s/c  gives 


and 


s _ n  -  a 

c  n|(H)  -  0T](H) 


(13) 


J  =  6(H)  -  Ti(H) 
6(H)  -  J  tj(h) 


(14) 


Recalling  that  H  is  determined  by  P^  for  a  given  6(Y)  curve  (given 
n,  y0'),  one  sees  that  equation  (14)  determines  I  as  a  function 
of  M.,  n,  a,  P-^,  and  y0’.  A  convenient  determination  of  I,'  of  course, 
can  only  be  made  if  the  function  tj(h)  in  addition  to  |(y)  has  been 
computed.-  The  function  t](h),  which  for  completeness  should  be  written 
as  tj(H;  M*,,  n,  yQ  ’),  is  somewhat  easier  to  compute  than  |(y)  since  it 
is  not  singular  at  Y  =  1.  Attention  is  called  to  the  fact  that  all  the 
above  integrals  with  limits  ranging  from  Y  =  0  to  Y  =  H,  as  in  equa¬ 
tion  (12),  for  example,  really  correspond  to  integration  over  both 
curved  surfaces,  first  from  Y  =  0  to  Y  =  1,  and  then  from  Y=1  to  Y  =  E 


With  the  position  of  maximum  thickness  determined  by  equation  (13), 
the  maximum  thickness  ratio  can  be  determined  in  terms  of  the  surface 
slope  at  the  leading  edge. 


'dY  \  _  2s  /dy\ 
dX  t  ylx )Q 


or,  from  equation  (ll). 


Jt 

c 


2 


yo1 

(dY/dX)0 


in  an  alternate  form 


►  (15) 

. 


t^ 

c 


2 

VV  in 


(15a) 


Calculation  of  pressure  drag  of  an  optimum  -profile.-  In  reference  4 
it  was  shown  that  for  linearized  supersonic  airfoil  theory  the  pressure 
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drag  coefficient  of  an  optimum  profile  was  a  simple  algebraic  function 
of  certain  quantities  such  as  H  and  s/c.  Since  these  quantities  are 
known  once  the  shape  of  the  optimum  profile  is  determined,  a  separate 
integration  is  not  required  in  order  to  calculate  the  pressure  drag. 
Fortunately,  a  similar  algebraic  relation  can  also  be  developed  for  the 
present  case.  In  so  doing,  integration  by  parts  is  employed  starting 
with  the  defining  equation  for  pressure  drag. 

|  cd  =  /*/2  Pdy  -  Pb  |  =  Pi§  ydP  -  Pb|  (16) 

> 

In  these  equations,  and  subsequent  ones,  the  integration  is  carried  out 
only  over  the  two  curved  portions  since  the  straight  midsection  can  con¬ 
tribute  no  drag.  From  equation  (3a), 


dP 


s(i/y)  d 


,  _y.2  |E-d 

y'  J  <V  \  AjJ 


- -x  [(t/2)n  -  yn)  a(j|) 


hence,  substituting  into  equation  (l6)  and  again  integrating  by  parts 
gives 


2  ca  2 


P-L  “  Pb  +  y. 


Ap 

W7 


1-1 


-\  f  l(t/2)n  -  (n  +  D/lax  - 


xfs+i  -  (n  +  dx 

The  first  bracketed  term  on  the  right  side  vanishes  by  virtue  of  equa¬ 
tion  (4).  The  remaining  integrals  can  be  simplified  by  noting  that 


I 


1  rcX- 

CJ0  TZW 


dx. 


c 


and  X,(t/2)n 


There  results 


ca  =  2yo,Z(^r)  t(n  +  1  -  a)T  -  1]  (17) 

This  equation  enables  the  pressure  drag  to  be  readily  calculated  if  the 
base  pressure  is  given,  since  P^,  determines  H  for  a  given  y0'  and 
Moo,  and  H  determines  T  in  accordance  with  equation  (14).  Thus, 
equation  (17)  involves  the  base  drag  implicitly,  but  not  explicitly. 

For  the  special  case  of  linearized  supersonic  flow, 


dP 

8y' 


=  2/V  Mro2  -  1 


and  the  above  equation  for  cd  can  be  shown  to  reduce  to  the  correspond 
ing  equation  for  pressure  drag  developed  in  reference  4. 
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*  Closed-Form  Solution  for  the  Special  Case  of  Given 

Cross-Section  Area  n  =  1,  cr  =  0 


When  the  cross-section  area  of  a  profile  is  prescribed  (n  =  1, 
o=0),  corresponding  to  a  given  wing  volume,  torsional  stiffness,  or 
torsional  strength  of  a  thin-skin  structure,  then  the  differential 
equation  (3)  can  be  integrated  immediately  with  respect  to  x  to  yield 
a  solution  in  closed  form  for  the  airfoil  shape.  There  results 


-\x  =  P  +  y*  sr-—  +  constant 
ay* 

The  constants  can  be  eliminated  by  evaluating  this  expression  at 
and  x  =  s  to  obtain 


(18) 


X  =  0 


* 


P  +  y^dp/dy1 )  -  P(p) 
po  +  V^P/dy^o  -  P(0) 


(19) 


Here  P(o)  is  the  pressure  coefficient  at  y»  =  0,  and  PD  is  the 
pressure  coefficient  at  x  =  0.  For  practical  purposes  P(o)  usually 
can  be  taken  as  zero,  although  strictly  speaking  it  should  be  regarded 
•as  a  small  quantity  compared  to  PQ  +  y0» (dP/dy*  )D.  The  parametric 
equations  for  Y  and  H  in  t.erms  of  y»  are  the  same  as  before,  only 
with  n  =  1. 


Y  _  1  _  y,g(5p/5y») 

y0,2(ap/ay«)0 

H  =  1  -  yi*  a(5p/9y»)1 
y0»2(8p/By*)0 


(20) 

(21) 


The  general  equation  for  the  base  pressure  coefficient  does  not  involve  n 
and  hence  is  the  same  as  before.  ’ 

.  ^  ■  Pl  +  (tp^  « 

The  constant  in  equation  (18)  can  be  evaluated  at  x  =  c  instead 
of  at  x  =  s.  Combining  such  an  evaluation  with  the  above  equation  for  Pb 
yields  the  alternate  expression  D 

x  _  1  P  +  y»(dp/&y»)  -  Pb 

c  p0  +  y0'(te/%r)0  -Pb 

which  involves  Pb  instead  of  P(o).  The  equations  for  t/c,  s/c,  and  cd 
can  be  derived  easily  from  the  preceding  equations.  Omitting  algebraic  ° 
details,  the  following  results  are  obtained: 
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By  evaluating  equation  (l8)  at  x  =  c,  and  combining  with  equations  (3b) 
and  (4), 

t  =  .  2yp,g (Sp/Sy* )0  (23) 

c  ~  P0  +  yo’^P/dyOo  -  pb 

by  evaluating  equation  (22)  at  x  =  s, 

s  _ _ P(0)  -  Pb  (24) 

c  P0  +  yo'^P/dy'Jo  -  pb 

and  by  substituting  n  =  1  and  a  =  0  into  equation  (17)* 

=a  =  2  y/  <*  -  1)  (=5) 

Solution  for  the  Special  Case  of  Given  Thickness  Ratio 

n  =  at,  a  =  finite 


In  reference  4  it  was  shown  that  the  limiting  values  n  =  00, 
a  =  finite,  represent  the  auxiliary  condition  of  a  given  airfoil  thick¬ 
ness  ratio.  The  mathematical  simplification  inherent  in  the  use  of 
approximate  theories  such  as  linearized  flow  enables  the  solution  for 
a  given  thickness  ratio  to  be  obtained  directly  by  passing  the  general 
solution  to  the  limit  as  n-»«.  For  shock-expansion  theory,  though, 
a  general  solution  in  closed  explicit  form  cannot  be  obtained,  and 
recourse  to  the  alternate  method  indicated  in  reference  4  is  required. 
This  alternate  method  deals  directly  with  the  appropriate  differential 
equation,  which,  for  the  case  of  given  airfoil  thickness,  becomes  simply 

y»2  =  constant  (26) 

ay' 

which  is  satisfied  by  a  profile  composed  of  any  number  of  straight  seg¬ 
ments.  As  shown  in  reference  4,  the  constant  in  the  above  equation  does 
not  change  over  the  entire  chord,  with  the  result  that  the  upper  half 
of  the  profile  forward  of  the  -trailing  edge  is  composed  of  two  straight 
lines,  one  extending  from  x  =  0  =  y  to  x  =  s,  y  =  -t/2,  and  the  other 
extending  from  x  =  s,  y  =  t/2  to  the  trailing  edge  x  =  c,  y  -.h/2. 

The  slope  is  discontinuous  at  the  point  where  y  =  t/2.  To  obtain  a 
solution  using  any  given  airfoil  theory,  it  is  necessary  to  satisfy  the 
differential  equation  (26),  the  end  condition  (4),  and  the  boundary 
condition  of  a  fixed  thickness  ratio. 

Equation  (17)  for  cd  becomes  indeterminate  as  n->  00  because 
7->0.  For  this  case,  however,  the  shape  is  known  and  the  pressure  drag 
be  determined  from  simple  physical  considerations: 

(cd) given  t/c  =  c  ^Po  "  M1  “  H)  PbH'1 


can 
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APPLICATION  OF  ANALYSIS 
Flow  at  Infinite  Mach  Number 


Prior  to  considering  shock-expansion  theory,  the  relatively  simple 
case  of  flow  at  infinite  Mach  number  over  slender  airfoils  with  small 
surface  curvature  will  be  considered.  For  such  conditions  the  pressure 
coefficient  on  a  surface  facing  upstream  is  proportional  to  the  square 
of  the  local  surface  slope,  that  is,  P  =  Cy'2.  Since  the  pressure  coef¬ 
ficient  is  zero  on  any  surface  facing  downstream  it  follows  from  physical 
considerations  that  H  =  1.  Equation  (4)  is  satisfied  by  requiring  that 
y ! *  =  0 .  By  consideration  of  the  differential  equation  (3a)  as  special¬ 
ized  to  the  present  case  it  follows  that 


1/6 


By  substitution  into  equation  (9),  and  employment  of  gamma  functions  to 
evaluate  the  integral  in  the  denominator,  there  results 


The  function  Y  versus  £(Y)  is  plotted  in  figure  3  for  n  =  .1,  2,  3,  and  ». 
The  infinite  value  of  n  corresponds  to  the  auxiliary  condition  of  a 
given  thickness  ratio,  and  the  optimum  profile  in  this  case  is  a  wedge, 

since  -I'i-)  ->1  and  £(Y)->Y  ns  n-»«.  It  is  seen  that  there  is 
n  \n/ 

little  difference  between  the  three  curves  for  finite  n. 


The  other  characteristic  function  needed  for  the  complete  determina¬ 
tion  of  an  optimum  profile  is  t}(h).  By  substituting  H  =  1  in  equa¬ 
tions  (12),  (13),  (14),  and  (17),  and  employing  equations  (5),  (10),  and 
(28),  the  following  expressions  are  obtained: 


n(i) 


2n 

2n  +  3 


—  =  (n  ~  cr)(2n  +  3) 
c  n(2n  +  3  -  2c) 


I 


3 

2n  +  3  -  2n 


(28) 

(29) 

(30) 


l _ 30 _ 

n(2n  +  3  -  20) 


(3D 


c 


Ik 
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2  /l\  n“a  2  f(3  +  2n  -2a) 


c  C 


_  S.  f- 

~  c 


n  -  a 


°3  -  iiCyo'  [2n  +  3  -  2o) 


/t\3  /  \3  f  (f^r  (n) 

=  w(§)  W  “TTT-r- 

■  ■  2  rV3  “/ 


(n  ~  g) 

(2n  +  3  -  2a) 


If  desired,  this  last  equation  for  the  pressure  drag  coefficient  can  be 
written  in  terms  of  I  instead  of  t/c,  inasmuch  as  I  is  related  to 
the  thickness  ratio  through  equation  (32).  It  should  be  noted  that  e 
above  equations  relate  in  closed  form  all  pertinent  properties  of  the 
optimum  profile  to  the  given  quantities  I,  n,  and  a.  Examples  of  op 
mum  profiles  determined  with  the  aid  of  these  equations  are  presented 

subsequently. 

Shock-Expansion  Theory 

When  the  oblique  shock-wave  and  Prandtl-Meyer  equations  are  combined 
to  calculate  the  pressure  on  an  airfoil  surface  in  supersonic  flow,  the 
resulting  equations  for  P  are  quite  involved.  The  appropriate  equa¬ 
tion  for  Bp/By’ ,  however,  may  be  obtained  by  starting  with  the  loca 
differential  relation 

1  5?  _ 2 _  (or,\ 


|pv2  as 


f M2  -  1 


This  point  relation  is  formally  the  same  as  the  corresponding  reiation 
applied  throughout  an  entire  flow  field  m  linearized  superso 
theory.  The  partial  derivative  Bp/B&  is  taken  with  M.  and  SQ 
constant.  Expressing  equation  (35)  in  terms  of  the  pressure  coefficient 
and  free-stream  conditions, 


or,  since 


BP  dP  2  , 

- —  =  - —  cos  o 

By’  BS 


there  results 


BP  0V2  2  cos2  S 
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*  This  equation  appears  fairly  simple,  but  is  not  in  a  form  which  can  be 
computed  readily  with  the  aid  of  existing  tables  where  quantities  such 
as  local  pressure  ratio  p/pt  and  total -pressure  ratio  across  a  shock 

*  wave  pt/pt  are  tabulated.  Thus,  a  more  convenient  form  for  calcula¬ 
ting  purposes  is 


From  this  equation  a  numerical  value  of  8p/dy'  can  be  determined  from 
tabulated  oblique-shock  and  expansion  characteristics  once  y0!,  Mw, 
and  the  local  slope  y*  are  specified. 

The  functions  £(y)  and  q(H)  have  been  calculated  for  shock- 
expansion  theory  by  substituting  equation  (36)  into  equations  (9)  and  (12), 
respectively,  artd  then  performing  the  indicated  integration  graphically 
by  the  method  outlined  in  appendix  A.  In  this  process  other  useful 
quantities  are  calculated  such  as  (dp/dy* )0  and  jn.  The  results  are 
presented  in  figures  4,  6,  and  7»  -  In  figure  4  the  quantity 

(dP/dy'  )o/(2A/m002  -  1),  which  is  equal  to  (dP/dy’  )0/[(Sp/Sy'  )o]Bo=0o> 

is  plotted  as  a  function  of  80  for  various  values  of  Mw.  Similarly, 
jn/(jn)6  =0o  is  plotted  in  figure  5*  It  is  to  be  noted 


Curves  of  Hn  versus  80  for  various  values  of  base  pressure  are  pre¬ 
sented  in  figure  6,  from  which  it  is  apparent  that  the  dimensionless 
trailing-edge  thickness  increases  if  either  the  airfoil  thickness  . 
increases  (S0  increases),  or  if  the  base  pressure  increases.  In  figure  7 
the  functions  |(Y)  and  q(H)  for  various  80,  and  n  are  Presented 
plotted  in  the  form  ,Y  versus  6,  and  H  versus  rj.  The  curves  of  Y 
versus  6  determine  the  shape  of  the  optimum  profile,  while  the  curves 
of  H  versus  tj  are  useful  in  determining  I,  s/c,  and  t/c.  Many  of  the 
curves  of  Y  versus  £  have  been  terminated  at  the  point  (indicated  by 
small  circle)  corresponding  to  zero  base  pressure. 


EXAMPLES  AND  DISCUSSION 


In  order  to  determine  an  optimum  profile  it  is  necessary,  of  course, 
to  know  the  base  pressure.  Experiments  have  shown  that  base  pressure  in 
two-dimensional  flow  depends  principally  on  the  Mach  number,  type  of 
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boundary- layer  flow,  and  the  boundary- layer  thickness  at  the  base.  * 

(See  reference  6.)  Average  experimental  values  are  shown  in  figure  8 
for  both  laminar  an4  turbulent  flow  plotted  as  a  function  of  the  parame¬ 
ters  proportional  to  the  ratio  of  boundary-layer  thickness  to  trailing- 
edge  thickness.  Shep-by-step  details  of  the  method  of  determining  an 
optimum  profile  by  combining  experimental  base  pressure  data  with  the 
curves  of  figures  4  to  7  are  given  in  appendix  B. 

In  figure  9  examples  of  optimum  profiles  determined  by  the  theory 
of  the  present  report  are  shown  together  with  corresponding  profiles 
determined  by  linearized  theory  (reference  4).  For  each  of  the  various 
auxiliary  conditions  the  particular  value  of  I  selected  for  these 
examples  is  equal  to  that  for  a  circular-arc  biconvex  airfoil  of  thick¬ 
ness  ratio  tca/c  =  0.06.  Since 

22n(n!)2(tca/2)n-g 
ca  (2n  +1)1 

(see  reference  4),  it  follows  that  with  c  =  1  the  optimum  profiles  in 
figure  9  correspond  to  the  value 

j  _  22n(n»)2(0«03)n~Cr 

(2n  +l)l  * 

The  auxiliary  condition  for  n  =  w  (fig.  9(d))  corresponds  to  a  given 

maximum  airfoil  thickness  of  0.06c.  As  indicated  in  figure  9  (and  also  -  ^ 

in  subsequent  figures  of  this  report),  the  base  pressure  for  M*  =  1.5 

and  Mj,,  =  3.0  corresponds  to  turbulent  boundary-layer  flow  at  a  Reynolds 

number  of  107.  Since  h  is  involved  in  the  abscissa  of  figure  8,  due 

allowance  is  made  for  the  variation  of  base  pressure  with  trailing-edge 

thickness.  Because  base  pressure  data  are  not  available  as  yet  for 

M,,  =  8,  a  constant  value  has  been  assumed  (pb  =  0.1  p^,)  which  is  believed 

to  be  reasonable  for  a  moderately  thick  trailing  edge  (h/c  0.05  or  more), 

but  probably  greatly  overestimates  the  base  drag  for  a  thin  trailing  edge 

(h/c  ~  0.01  or  less).  For  =  00,  it  is  not  necessary  to  know  the  base 

pressure  since  the  optimum  profile  at  this  limit  is  independent  of  Pb» 

In  figure  10  examples  are _ shown  for  various  values  of  I  with 
M,*  =  3,  n  =  3,  and  a  =  0.  Instead  of  specifying  the  value  of  I  in  each 
case,  the  thickness  ratio  of  a  structurally  equivalent  circular-arc 
biconvex  airfoil  is  specified,  as  the  significance  of  this  latter  value 
is  easier  to  visualize.  As  would  be  expected,  there  is  no  appreciable 
difference  between  the  profiles  determined  by  linearized  and  shock- 
expansion  theory  when  the  value  of  I  is  small  (tCa/c  =  0.02,  for 
example),  although  differences  are  evident  for  larger  values  of  I 
(tca/c  =  0.04  and  larger). 

For  each  of  the  examples  shown  in  figures  9  and  10,  it  is  to  be 
noted  that  in  comparison  to  the  profile  determined  by  linearized  theory, 
the  corresponding  optimum  profile  determined  by  shock-expansion  theory 
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has  a  smaller  slope  over  the  portion  of  surface  facing  upstream,  a 
position  of  maximum  thickness  farther  aft,  and  a  greater  slope  over  the 
portion  of  surface  facing  downstream.  This  difference  which  increases 
with  increasing  Mach  number  is  to  be  expected,  as  indicated  in  refer¬ 
ence  4,  because  the  linearized  theory  overestimates  the  suction  forces 
and  underestimates  the  positive  pressure  forces. 

As  regards  drag  coefficient,  it  is  evident  that  drag  calculations 
based  on  linearized  theory  cannot  be  used  at  hypersonic  Mach  numbers 
since  the  computed  coefficient  approaches  zero  as  the  Mach  number 
increases.  Also,  it  is  to  be  remembered  that  linearized  theory  is  con¬ 
siderably  less  accurate  in  predicting  the  drag  of  blunt-trailing-edge 
airfoils  than  of  sharp-trailing-edge  airfoils,  since  the  Busemann  second- 
order  terms  for  the  upstream  and  downstream  surfaces  do  not  cancel  as 
they  do  when  the  trailing  edge  is  sharp. 

If  shock-expansion  theory  is  used  to  calculate  the  drag  of  the  pro¬ 
file  determined  by  linearized  theory,  the  resulting  value  is  only  slightly 
greater  than  the  drag  of  the  same  profile  determined  by  shock-expansion 
theory.  In  order  to  put  this  idea  on  a  more  firm  quantitative  basis, 
the  case  of  infinite  Mach  number  can  be  considered,  as  the  differences 
between  profiles  determined  by  linear  and  nonlinear  theory  are  the 
greatest  at  this  limit.  (See  examples  in  fig.  9.)  By  use  of  the  expres¬ 
sion  P  =  Cy*2  to  calculate  the  drag  of  the  profile  determined  by  linear 
theory,  and  dividing  by  the  drag  calculated  correspondingly  for  the  pro¬ 
file  determined  by  nonlinear  theory,  the  following  expression  results: 

1.08  for  n=l,q=0 
1.06  for  n=2,a=0 
1.07  for  n=2,cr=l 
1.04  for  n=3,a=0 
1.05  for  n=3,cr=l 
1.00  for  n=co,a  finite 

(38) 

It  is  seen  that  at  infinite  Mach  number  the  actual  pressure  drag  of  an 
optimum  profile  whose  shape  is  determined  by  linearized  theory  does  not 
exceed  the  drag  of  the  true  optimum  profile  by  more  than  about  8  percent. 

If  consideration  is  given  to  the  consistent  differences  noted  earlier 
between  the  shapes  of  optimum  profiles  determined  by  linear  and  nonlinear 
theory,  it  is  evident  that  linearized  theory  can  be  used  with  good 
accuracy  to  determine  the  optimum  profile  at  any  supersonic  Mach  number 
up  to  infinity.2  Even  without  considering  the  consistent  difference  noted 
above,  the  profile  determined  by  linearized  theory  is  sufficiently  accurate 
for  most  engineering  purposes.  Under  less  general  conditions  a  similar 

2As  indicated  in  reference  4,  however,  the  linearized  theory  does  not 
yield  a  reasonably  accurate  profile  at  the  low  supersonic  Mach  numbers 
near  or  below  shock  detachment.' 
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result  also  has  been  found  in  the  recent  investigation  of  Klunker  and 
Harder  (reference  7)  which  appeared  while  the  present  report  was  being 
prepared.  The  shape  of  some  of  the  optimum  profiles  determined  in 
reference  7,  however,  doe 3  not  agree  with  the  shape  of  analogous  profiles 
in  this  report.  For  example,  it  is  indicated  in  reference  7  that  the 
profile  of  least  drag  for  t/c  =  0.06  has  a  sharp  trailing  edge  at  all 
Mach  numbers  below  about  6,  whereas  the  corresponding  profiles  shown  in 
figure  9(d)  indicate  appreciable  trailing-edge  thickness  even  at  Mach 
numbers  of  1.5  and  3.  This  discrepancy  is  attributed  to  the  arbitrary 
base  pressure  curve  assumed  in  reference  7  which  does  not  correspond  to 
measured  data  for  thin  trailing  edges. 

From  an  engineering  viewpoint  it  is  desirable  to  know  how  much  lower 
the  drag  of  an  optimum  profile  is  than  that  of  a  sharp-trailing-edge 
profile,  and  also  how  much  the  optimum  profile  can  be  altered  without 
significantly  increasing  the  drag.  In  order  to  provide  a  basis  of  com¬ 
parison,  the  zero-lift  pressure  drag  of  a  family  of  sharp-trailing-edge 
circular-arc  biconvex  airfoils  of  various  thickness  ratios  has  been 
calculated  by  shock-expansion  theory  for  the  Mach  number  range  between 
1.5  and  8.  The  results  are  shown  in  figure  11.  Thus,  for  any  profile 
the  drag  of  a  structurally  equivalent  (same  value  of  i)  circular-arc 
biconvex  profile  can  be  determined  readily  from  the  curves  in  figure  11 
by  simply  calculating  tCa  from  the  equation 

I  -  lea  -  22n(n!)2(tea/2)n-°/(2n  +1)! 

Computations  of  drag  have  been  made  for  a  family  of  "semioptimum" 
profiles  having  arbitrarily  selected  values  of  trailing-edge  bluntness  H, 
a  shape  forward  of  the  trailing  edge  that  yields  minimum  foredrag  for  each 
particular  H,  and  the  same  value  of  I  as  a  circular-arc  biconvex  pro¬ 
file  of  thickness  ratio  tea.  These  calculations  have  been  carried  out 
for  tqa/c  =  0.02,  0.04,  0.06,  0.08,  and  0.10  at  Mach  numbers  of  1.5*  3> 
and  8,  and  for  various  combinations  of  n  and  cr.  As  in  previous 
examples,  the  base  drag  in  each  case  was  determined  from  the  curves  of 
figure  8  for  turbulent-boundary- layer  flow  at  He  = :  107.  The  results  are 
shown  in  figure  12  plotted  in  the  form  of  a  drag  ratio  versus  H.  Each 
curve  corresponds  to  a  constant  value  of  I,  and  is  identified  by  the 
thickness  ratio  (tca/c)  of  a  circular-arc  biconvex  profile  having  the 
same  value  for  I.  In  order  to  maintain  a  constant  value  of  I,  the 
actual  thickness  ratios  (t/fc)  of  the  semioptimum  profiles  change  somewhat 
as  H  varies  between  0  and  1  (the  ratio  t/tca  lies  between  about  0.90 
and  1.05  for  the  case  n  =  1,  a  =  0,  between  about  0.71  and  0.84  for 
n  =  2,  o  =  1,  and  between  about  0.97  and  1.08  for  n  =  3>  cr  =  0).  For 
each  curve  in  figure  12  the  semioptimum  profile  having  the  minimum  drag 
coincides  with  the  optimum  profile  determined  from  the  curves  of  fig¬ 
ures  5  to  8.  The  ordinate  of  each  minimum  point  indicates  the  relative 
drag  of  the  optimum  compared  to  a  structurally  equivalent  circular-arc 
biconvex  profile,  while  the  rise  on  each  side  of  the  minimum  indicates 
the  drag  penalty  resulting  from  the  use  of  too  much  or  too  little 
trailing-edge  thickness.  It  may  be  noted  that  some  of  the  curves  do  not 
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cover  the  complete  range  of  values  of  H.  In  all  such  cases,  however, 
sufficient  calculations  were  made  so  that  the  minimum  point  was  bracketed. 

Perhaps  the  most  significant  feature  apparent  from  the  curves  of 
figure  12  is  the  large  increase  in  drag  that  results  for  thin  airfoils 
when  a  trailing-edge  thickness  considerably  greater  than  the  optimum  is 
employed  at  moderate  supersonic  Mach  numbers.  The  pressure  drag  of  full- 
blunt  profiles  at  Mach  numbers  of  1.5  and  3  for  the  case  tca/c  =  0.02  is 
several  times  the  pressure  drag  of  the  optimum,  the  exact  factor  varying 
between  about  2.8  and  3.9  depending  on  M^  n,  and  cr.  On  the  other 
hand,  for  thicker  profiles  (tca/c  aj  0.08)  the  pressure  drag  is  much  less 
sensitive  to  variations  in  trailing-edge  thickness  about  the  optimum, 
and  the  use  of  a  full-blunt  profile  instead  of  the  optimum  would  result' 
in  a  much  smaller  percent  drag  penalty.  Also,  for  a  given  value  of  tca/c 
it  is  evident  from  figure  12  that  the  pressure  drag  becomes  less  sensitive 
to  variations  in  trailing-edge  thickness  from  the  optimum  as  the  Mach  num¬ 
ber  is  increased  into  the  hypersonic  regime. 

As  would  be  expected,  the  semioptimum  sharp-trailing-edge  profiles 
(H  =  0  in  fig.  12)  have  somewhat  less  pressure  drag  than  a  structurally 
equivalent  circular-arc  biconvex  profile.  The  observed  difference  in 
drag  for  thin  airfoils  at  moderate  Mach  numbers  is  negligible  for  the 
case  of  n  =  1,  0  =  0,  since  the  optimum  sharp-trailing-edge  profile  for 
these  conditions  is  very  close  to  a  circular-arc  biconvex  profile  (if 
linearized  theory  were  employed  the  optimum  sharp-trailing-edge  profile 
for  n  =  1,  0=0  would  be  a  circular-arc  biconvex  profile).  The 
corresponding  difference  in  drag  for  the  case  n  =  2,  0=1,  however,  is 
significant  since  the  optimum  profile  in  this  case  has  a  midsection  of 
constant  thickness,  and  hence  is  of  considerably  different  shape,  as 
well  as  being  considerably  thinner  than  a  structurally  equivalent  circular- 
arc  biconvex  profile. 

Each  of  the  curves  for  Mb*  »  1.5  and  3  (figs.  12(a)  and  12(b))  show 
a  minimum  at  some  finite  value  of  H,  but  the  curves  for  tca/c  =  0.02  at 
M  =  8  (fig.  12(c))  do  not.  If  all  other  parameters  were  constant,  this 
trend  would  not  be  expected  inasmuch  as  the  optimum  trailing-edge  blunt¬ 
ness  for  a  given  thickness  ratio  generally  increases  as  the  Mach  number 
is  increased  into  the  hypersonic  range.  (See  reference  4.)  The  unexpected 
trend  is  observed  in  the  present  examples  because  the  variation  of  base 
pressure  with  trailing-edge  thickness  is  considered  at  M  =  1.5  and  3 
where  experimental  measurements  are  available,  but  it  is  not  considered 
at  Mo#  =  8  where,  in  the  absence  of  experimental  data,  a  constant  base 
pressure  was  arbitrarily  assumed  (one-tenth  of  the  free-stream  pressure, 
irrespective  of  trailing-edge  thickness).  It  is  expected  that  if  base 
pressure  measurements  were  made  at  M*,  =  8,  they  would  show  a  dependence 
on  trailing-edge  thickness  just  as  at  the  lower  Mach  numbers.  Conse¬ 
quently,  it  is  believed  that  the  actual  curves  for  the  thinner  airfoils 
at  M,,  =  8  will  be  greatly  different  than  shown  in  figure  12(c),  although 
the  curves  for  the  thicker  airfoils  are  not  expected  to  be  significantly 
different.  If  a  constant  base  pressure  corresponding  to  measured  values 
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on  thick  trailing  edges  were  assumed  at  M^,  =  1.5  and  3?  the  curves 
for  tca/c  =  0.02  would  rise  starting  from  H  =  0  even  more  steeply 
than  at  M^  =  8.  This  illustrates  the  necessity  of  considering  the 
dependence  of  base  pressure  on  trailing-edge  thickness  in  such  an 
analysis.3 

The  drag  of  the  various  optimum  profiles  (minimum  points  in  fig.  12) 
is  seen  to  be  less  than  the  drag  of  a  structurally  equivalent  circular- 
arc  biconvex  profile  by  amounts  varying  between  about  1  percent  and 
53  percent,  depending  on  the  values  of  n,  0,  M^,  and  tca/c.  The  largest 
drag  reduction  occurs  for  the  case  of  n  =  2,  a  =  1,  just  as  indicated 
by  the  linearized  theory  of  reference  4.  Likewise,  the  drag  reduction 
generally  increases  as  the  Mach  number  or  the  thickness  ratio  is 
increased,  as  predicted  by  linearized  theory. 


CONCLUSIONS 


1.  For  a  given  Mach  number  and  structural  requirement,  the  shape 
and  drag  of  the  profile  having  the  least  possible  pressure  drag  at  zero 
lift,  as  computed  from  shock-expansion  theory,  can  be  determined  readily 
provided  the  base  pressure  is  known  from  experiments  and  provided  curves 
of  certain  characteristic  functions  are  available.  (These  functions 
are  |(Y,  6d,  n),  q(H,  60,  n),  H(60,  Pb/Pco,  n),  jn(&o)>  and  8p/8y»  (bD)  > 
examples  of  which  are  presented  in  figures  4  to  7.) 

2.  A  comparison  of  profiles  determined  by  shock-expansion  theory 
and  linearized  theory  indicates  that  the  linearized  theory  may  be  used 
with  reasonable  accuracy  at  Mach  numbers  up  to  infinity  to  determine  the 
shape  of  the  optimum  profile,  although  it  can  be  used  only  at  moderate 
supersonic  Mach  numbers  to  determine  the  drag. 

3.  Considerable  deviations  in  profile  shape  from  the  theoretical 
optimum  can  be  made  without  increasing  the  drag  excessively  provided  the 
Mach  number  is  high,  or  the  airfoil  thickness  ratio  is  relatively  large. 
Large  drag  penalties  result,  hpwever,  if  a  trailing-edge  thickness 
appreciably  greater  than  the  optimum  is  employed  on  a  thin  airfoil  at 
moderate  supersonic  Mach  numbers. 

4.  It  is  necessary  to  consider  the  experimentally  observed  depend¬ 
ence  of  base  pressure  on  trailing-edge  thickness  when  calculating  the 
opt imum -prof ile  shape  and  drag  of  a  thin  airfoil. 


Ames  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 

_  Moffett  Field,  Calif.,  June  25,  1952  _  .  _ _ 

3lf  the  boundary. layer  were  laminar,  the  effect  on  the  optimum  shape  of 
the  dependence  of  base  pressure  on  trailing-edge  thickness  would  be 
even  greater. _ _ _ _ _ 
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APPENDIX  A 


METHOD  OF  CALCULATING  |(Y) 


The  singularity  at  Y  =  1  of  the  integrands  in  equation  (9)  causes 
difficulty  when  directly  evaluating  such  integrals  numerically  or  graphi¬ 
cally.  This  difficulty,  however,  can  be  circumvented  by  transforming 
from  Y  as  the  integration  variable  to  a  new  function  fn(Y)  defined  as 
follows : 


r*(Y)  7r 


With  this  transformation  the  equation  for  l(Y)  becomes 

rfn«/>r 


!(Y)  = 


/n  yli 

_  X  .  J  ay  n 


^  /cjp" 


In  this  equation  the  integrands  and  the  ranges  of  integration  are  all 
finite.  The  constant  kn  is  the  same  quantity  as  that  used  in  refer¬ 
ence  4,  namely. 


kn  =  fn(l)  = 


for  n  =  1 


for  n  =  2 


1.4023  for  n  *  3 


for  n  =  00 


The  functions  fn(Y),  apart  from  an  additive  constant,  are  likewise  the 
same  characteristic  functions  as  appeared  throughout  the  analysis  of 
reference  4  when  linearized  supersonic  airfoil  theory  was  employed  for  P. 
Thus, 

f 2(1  -  Vl  -  Yn)  for  n  =  1  1 


fn(Y)  =< 


1.4023  -  3"4F(k,<P) 


for  n  =  2 


for  n 


for  n  =  00 
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where  F(k,9)  is  the  incomplete  elliptic  integral  of  the  first  kind  of 
modulus  k  =  sin  75°  =  0.9659*  and  amplitude 


Cp  = 


cos  1 


C/T  -  1  +  Y) 

C/T  +  1  -  y) 


With  the  transformation  to  fn,  the  integrals  in  equation  (A2)  are 
evaluated  by  first  selecting  a  number  of  values  of  y‘  ranging  from  y0* 
to  large  negative  values.  For  each  y*  the  ordinate  Y  is  computed 
from  equation  (6),  f^  from  equation  (A4) *  and  cSP /6y '  from  the  particu¬ 
lar  airfoil  theory.  A  plot  is  then  made  of  dP/dy'  versus  fn  in  order 
to  evaluate  the  integrals  determining  i(Y). 
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APPENDIX  B 

DETAILS  OF  METHOD  OF  DETERMINING  THE  SHAPE  AND  DRAG  OF 
AN  OPTIMUM  PROFILE  BY  SHOCK-EXPANSION  METHOD 


In  the  shock-expansion  equations  the  leading-edge  deflection 
angle  S0  is  a  more  convenient  parameter  to  use  than  the  given  value 
of  I,  hence  the  steps  outlined  below  involve  an  iterative  procedure. 

(1)  Assume  values  of  80  and  p^/p^ 

(2)  Read  value  of  Hn  from  figure  6  and  compute  H 

(3)  Read  values  of  6(H)  and  t}(H)  from  figure  7;  compute  T  from 

equation  (l4)  and  s/c  from  equation  (13) 

(A)  Read  value  of  (dp/dy’)0  from  figure  4,  jn/(jn)&0  =  o° 

from  figure  5;  compute  (jn)g  _  Qo  from  equation  (37), 
t/c  from  equation  (15a) 

(5)  Compute  I  =  T  (t/2)n“0 

By  comparison  of  the  computed  value  of  I  with  the  given  value,  a  new 
value  of  B0  can  be  estimated.  Also,  from  the  computed  value  of  h/c, 
the  experimental  base  pressure  curves  in  figure  8  yield  a  new  value  of 
Pb/p^.  By  repetition  of  the  above  steps  until  the  final  computed  value 
of  I  is  equal  to  the  given  value,  and  the  final  computed  value  of  h/c 
corresponds  to  the  final  base  pressure  assumed,  all  characteristics 
(So,  t/c,  s/c,  H,  S(Y),  (dP/dy»)  o>  jn>  an&  I)  of  the  optimum  profile  are 
determined.  The  pressure  drag  is  then  calculated  from  equation  (17). 
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y‘  -  tan  8 
yoston80 

y}=  tan  8,  A  B 


Figure  /.-  Sketch  illustrating  various  symbols  used  in  the  analysis. 


Figure  2.  -  Sketch  of  a  typical  curve  for  determining  optimum-airfoil 

profiles. 
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Figure  4.- Curves  based  on  shock-expansion  theory  showing  (—,)  /(2/JM&-1) 
as  a  function  of  leading-edge  deflection  angle  for  various  Mach  numbers. 


h/(h)s,.0’  J> /(*•)». 
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S0,  deg. 


Figure  5. — Curves  based  on  shock-expansion  theory  showing  jn/(jnh0*o' 
as  a  function  of  leading -edge  deflection  angle  for  various  Mach 
numbers. 
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S0,  deg 

(a)  Mco-1.5,  2,3 

Figure  6.  Curves  based  on  shock-expansion  theory  showing  dimension¬ 
less  trailing— edge  thickness  os  a  function  of  base  pressure  and 
leading-edge  deflection  angle. 
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c/[h(Re)'/z] 


(a)  Laminar  flow 


Figure  8.- Average  experimental  values  of  base  pressure  from  reference  6. 
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shock  expansion  theory 
linearized  theory 
approximate  theory  for  M^OO 


x/c 

(a)  n=/t  cr=0 

Figure  9 Examples  of  optimum  profiles }  auxiliary  integral  equal  to  that  of 
a  circular -arc  biconvex  airfoil  of  thickness  ratio  0.06. 


Shock  -expansion  theory 
Linearized  theory 
Approximate  theory  for  M«,=oo 


(b)  n  =  2,  <r  =1 
Figure  9.-  Continued. 
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Shock -expansion  theory 
Linearized  theory 


x/c 

(d)  n=oo,  <r  =  finite 
Figure  9-  Concluded. 


oK  n|s; 
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Shock-expansion  theory 
linearized  theory 


Figure  10-  Examples  of  optimum  profiles  for  values  of  auxiliary  integral  equal 
to  that  of  a  circular- arc  biconvex  airfoil  of  maximum  thickness  tea-,  n=3 , 
<r-0-,  Mce  =  3}  base  pressure  corresponding  to  turbulent  flow  at  Re- 1 07. 
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Figure  1 1  -Pressure  drag  of  circular-arc  biconvex  airfoils  determined  by  shock- 

expansion  theory. 


(a) M0 ‘1.5,  pb  corresponding  to  turbulent  flow  at  Re- 1 07 

Figure  12.- Effect  of  trailing- edge  bluntness  on  the  relative  pressure  drag 

of  semioptimum  profiles  compared  to  a  structurally  equivalent  circular- 
arc  biconvex  profile. 


Cd/(cd)ca  cd/(°dh 
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